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Abstract 


In the first part of the article, a new interesting system of difference equations is 
introduced. It is developed for re-rating purposes in general insurance. A nonlinear 
transformation ip of a, d-dimensional {d > 2) Euclidean space is introduced that 
enables us to express the system in the form := t = 0, 1, 2, .... Under 

typical actuarial assumptions, existence of solutions of that system is proven by 
means of Brouwer’s fixed point theorem in normed spaces. In addition, conditions 
that guarantee uniqueness of a solution are given. 

The second, smaller part of the article is about Leslie-Gower’s system of d > 2 
difference equations. We focus on the system that satisfies conditions consistent 
with weak inter-specific competition. We prove existence and uniqueness of the 
equilibrium of the model under surprisingly simple and very general conditions. 

Even though the two parts of this article have applications in two different 
sciences, they are connected with similar mathematics, in particular by our use of 
Brouwer’s Fixed point Theorem. 

Key words: Brouwer fixed point theorem. Loss ratio method, Leslie-Gower model, 
Beverton-Holt equations 
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1 Introduction and Preliminaries 

1.1 Introduction: To understand the purpose of the research let us first take a brief 
look into existing rating algorithms published in actuarial papers, even though it is not 
necessary for understanding mathematics of this article. 

The variety of rating algorithms can be found in [nms]. For example, Algorithm 1, 
in is derived from Bailey’s minimum bias conditions which in the case of two vector 
variables x and y with m and n dimensions (categories), respectively can be written as 
the following system of equations: 


n 



m 



Here, are observed loss costs and Wij are earned exposures to risk (weights) for the 
cells (*, j). The system has m + n equations with the same number of unknowns Xi and 
2 /j called relativities. 
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The algorithm derived from Bailey’s set of minimum bias conditions is the first and 
most important rating algorithm presented in It was first introduced in py but 

with four rather than two classification variables. 

Generally speaking for all algorithms introduced in papers mw, the authors started 
from a vector equations of the form 


</>(f)=0, (1.1) 

where the function (j) : D^f, integer d G N denotes generic dimension of an Euclidean 

space and C is domain of (j). Then they derived the corresponding vector equations 
of the form 

f = <p(f) (1.2) 

where ip : Dy, —>■ R‘^ and Dy, C R'^ is domain of ip. In the next step the corresponding 
algorithms are introduced by the formula 

f*+i :=(^(f‘),t = 0, 1, 2, .... (1.3) 

Unfortunately, in those actuarial papers the issues of necessary and sufficient conditions 

for convergence of those algorithms were not addressed. Another, actuarial rather than 
mathematical issue with algorithms from [D[SJ[5] is that neither follows standard and 
intuitive rating logic of adjusting rates to given experience. (The standard rating method 
is called loss ratio method. To learn about LR method, one might read e.g. [2114].) 

In this article, we will address both of the above mentioned shortcomings of the 
existing actuarial algorithms. We will derive an iteration formula of the form m 
directly from the standard loss ratio method. That was not possible to do earlier because 
the standard re-rating process was not mathematically formalized until relatively recently 
in [2]. The mathematical formalization enabled us to reduce number of inputs in the 
rating process; we do not need premium input any more. We will repeat the main 
features of the mathematically formalized rating model in Section [21 for convenience of 
the reader. 

Our algorithm is defined by system (EH) below. In Theorem 13.21 we prove that under 
typical assumptions used in rate-making area, system (13.11) has at least one solution. The 
proof is based on Brouwer’s fixed-point theorem. Sufficient conditions for convergence 
of algorithm EH towards an uniquely determined point are very likely satisfied under 
typical actuarial circumstances listed in Theorem l3.2l However, we have to add additional 
conditions in order to prove uniqueness of a fixed point of the function ip, cf. Theorems 
13.41 and and Corollaries 13.61 and 13.71 

In Section m we study a system of d > 2 difference equations that model competition 
of d species in an ecological system. This is a well-known Leslie-Gower (LG) model. Here 
we focus on LG model that satisfies certain conditions consistent with weak inter-specific 
competition. Similarly to the first part of this article, we use Brouwer’s fixed point theo¬ 
rem to prove existence of the solution (equilibrium) in the system. Then we give simple 
general conditions that guarantee uniqueness of the solution. We also prove that a linear 
algebraic system with the same solution can be joined to a given LG system of difference 
equations. Therefore, solving such LG systems simplifies down to solving linear algebraic 
systems. 

1.2 Preliminaries: We will here recall some basic concepts of functional analysis 
regarding Euclidean space R‘^ for the convenience of the reader, in particular for rating 
actuaries. 

Eor f = (/i,..., /d) e the functions defined by 

Ilf lip := (/f + ... + /D^^P,P>l, 
l|f|li:= l/i| + -.- + l/.|, 
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l|f|loo :=max{|/i|,...,|/d|}, 

are norms in Note that for p = 2 we get the Euclidean norm. 
It is known that the following relation holds 


l|flL<llf|lp<llf|li- 

It is also known that all norms defined in R‘^ are equivalent. It means that the convergence 
of the sequence {f*} C R‘^ towards f in one norm implies the convergence towards the 
same vector in other norms. Hence, we can deal with the most convenient norm for us 
in a particular situation. 

For given norms in i?" and i?™, the norm of a linear mapping M : i?" —>■ i?™ (which 
always can be represented by a matrix) is defined by 


l|M|| 


l|Mf|| 

l|f|| 


As all norms in R‘^ are equivalent and as matrix norms ||AI||j^ and ||M||^ have rather 
simple representations, we will focus on ||.||j^ and ||.||oo in this paper. 

Recall that for a function ip : Dip R‘^ {Dp C R'^) differentiable at x e Dp, 
Jacobi-matrix at x is 


Jp (x) — 


d{ipi,...,pd) 


/ 9yi (x 

dx\ 


dxd 


\ 


dx\ 


dxd 


d{xi, ...,Xd) 

Then, from well known formulae for matrix norms ||.||j^ and we get 

}, 

\Jv WIL = . ,{ -f ... -b }■ 


II'^V Will = .max { 

j = l,...,d 




WiW 

-b . 

. -b 

WdW 

dxj 



dx, 

WzW 



w*w 

9a; 1 



dxd 


Recall that by definition U C is a convex set if and only if 


(1.4) 

(1.5) 


(f, g G U, t G [0, 1]) — >■ tf -b (1 — t) g G U. 


2 Rating Model 

In the rate-making area all presentations are performed by means of two, three or, at 
most, four risk classification variables (or risk factors). We use three factors rather then 
two because generalizations to N > 3 factors are then straightforward. The three risk 
factors are presented here by vectors 

X = {xi, X2,... , Xm), y = (yi,?/2,--- , 2 /n), z = {zi, Z2, . ■ . , Zp) , 

where, by nature of the problem, 

x^ > 0, yj > 0, Zk > 0, i = I,... ,m-, j = 1,... ,n-, k = 1,... ,p. (2.1) 

Let us denote by D° C set of points that satisfy (12.ip . 

Risk factors determine the risk space, i.e. the set of all risk cells {i,j,k), where the 
triplet of risks quantified as {xi,yj,Zk) is assigned to every cell {i,j,k). If we keep one 
index fixed, we get a plane (slice) in the risk space. The base cell is usually the cell with 
the largest exposure to the risk (weight), so that it has maximal statistical credibility. 
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Without loss of generality we assume that it is the cell (1,1,1) and that the corresponding 
rate, the base rate = rm. The model is multiplicative. That means that the rate in the 
cell characterized by triplet of risks {xi,yj,Zk) is calculated by 

'^ijk — 'k’liiXiyjZ}^^, ( 2 . 2 ) 


Hence, we assume xi = yi = zi = 1, without loss of generality. 

The following algebraic representations of rating variables were introduced in the 
actuarial article [ 5 ]. The generic risk factors x, y, z were named there class, territory, 
industry, respectively. The coordinates of the risk factors were called, as usually, relativ¬ 
ities, or dijferentials. 

Let L — {lijk)^,^n,,,pilijk > 0, denote the three-dimensional array of fully developed 
and trended losses. We say that lijk is expected dollar loss or simply projected loss in the 
cell (i,j, k). Note that projections Ujk are calculated based on previous experience. The 
total expected (or projected) loss is, obviously, 

m n p 

2—1 j—1 k—1 

Three-dimensional array E = {^ijk)^ ft Oj consists of earned exposures to 
risk (weights) by cells. In actuarial terms, is number of units of insurance sold in the 
cell {i,j,k). Note that arrays L and E represent information about the past experience 
of the business, with rates calculated by (1^ . Note also that according to the nature of 
the problem, because exposure to risk is a necessary condition for loss, it holds 

lijk ^ 0 t eijk ^0. (2.3) 

The converse implication does not hold because it is possible to have positive exposure 
to risk in a cell {i,j, k), i.e. > 0 and that loss will not occur, i.e. lijk = 0. 

The purpose of the rating model is to adjust rates according to the information given 
by L and E. In order to do that we first adjust risk factors according to the information 
given by L and E. The adjusted risk factors are denoted by: 

X = (ii, X2,... , Xm) , y = {yi, y2,-..,yn), Z = {ZiZ2... Zp) 

and called indicated factors. Then indicated rates are denoted by and they are 
calculated by 

fijk = riiix^yjZk- (2.4) 

Indicated factors are calculated by formulas (12.8|) below and indicated base rate by 
formula (1^ . To understand those formulas we need to introduce some concepts. 

By keeping fixed one index at a time and summing loss amounts of the corresponding 
slice we obtain the vector of losses for factor x, denoted P, by 

n p 

^^=EE^bfc> i = l,2,...,m. (2.5) 

j=i k=i 

Vectors of losses for factor y, denoted P, and factor z, denoted E, are introduced 
similarly. Obviously, 

m n p 

S=1 S=1 s — 1 

Factor x adjusted exposures to risk are defined by 

n p 

^^=EE ^ijkyjZk^ i — 1, 2, ..., rn. (2.6) 

j=i k=i 
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Adjusted exposures and are introduced similarly. 

Factor x loss costs adjusted for heterogeneity, or simply factor x adjusted loss costs 
are defined by: 


T X ._ I 

i — ’ 


1 = 1,2, ... m. 


(2.7) 


and Lj and are introduced similarly. Then the following formulas for indicated factors 
were derived in Section 3.3 of [5]. 


Xi = 

Lf 

Lf’ 

i = 

1, 2, .. 

., m. 

Vj = 

LV 

j = 

1, 2, . 

.., n. 

Zk = 

LI 

LV 

k = 

1, 2, . 

.., p. 


( 2 . 8 ) 


Those formulas calculate the same indicated factors x, y and z as the standard loss ratio 
method, but with only losses and exposures as inputs, which is very helpful to rating 
actuaries. Namely, previously also premiums were needed for calculations of indicated 
factors and rates, see [3]. 

The following formula for indicated base rate has also been proven in [2] 


f^iii = 


L 1 


PER ^ " P _ , 

XrVjZkCijk 

i—1j—1k—1 


(2.9) 


where constant PLR stands for “permissible loss ratio”. The formulas (j2.4l) and (j2.9l) 
mean that indicated rates are calculated in terms of indicated factors, where the factors 
are obtained by means of simple formulas (12.81) that depend only on exposures and losses. 


3 Rating Algorithm 

3.1 System of difference equations: We now consider the formulas (12.81) to be merely 
the first iteration of the rating algorithm. For the second and all other iterations, we will 
use the same experience of losses L = and exposures E= {eijk)^,j,nxp- 

We could denote current risk factors by x°, y° and z° and the relativities after t*^ 
iteration by x*, y* and z*. However, we will avoid that notation in order to simplify this 
presentation. We will instead consider that formulas (12.81) calculate the next iteration, 
taking the previous one as input. Hence, in the following formulas, Xi, yj and Zk are 
values of t*^ iteration, while Xi, jjj and Zk are values of the next, {t + 1)*^ iteration. 

From (|2.6I) . (12.71) and (12.81) we obtain 


Vj = 


n p 

X) X) GljkVjZk 
If i=l k=l 

j^x n p ’ 

^ijkVj^k 

j^l k^l 
m p 

^ilk^i^k 

j i—\k—1 

jV m p ’ 

'EE ^ijk^i^k 

i—1k—1 


(3.1) 
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Zk 


,2 S S 

i-k _ 

j^z m n 

^ X/ S SijkXiUj 
i=lj=l 


In the next iteration, we substitute Xi, yj and Zk for Xi, yj and Zfe, into the formulas 
m, while values of the coefficients If, l^, If and Cijk, i = 1, m, j = 1, m, 
k = 1, ..., p, remain the same. The iteration process can be repeated until the iterations 
become sufficiently close, which we denote as 


x^^x^, i = l, yj ^ yj, j = 1, ..., n, Zk ^ Zk, k = 1, ..., p, 

where we define when the two iterations are close enough. If the iterative process con¬ 
verges, say to Xi, yj, Zk, then we can calculate the final indicated rates fijk by the 
formulas (1^ . (EH)' and (1^ . 


3.2 Existence of solutions: In this section we will prove that rating algorithm 
EH converges to a uniquely determined solution in many typical actuarial situations. 
Because all norms in Euclidean space are equivalent we will talk about convergence 
without mentioning any particular norm in the sequel. Let us first recall Brouwer’s 
fixed-point theorem which we restate as the following lemma. 

Lemma 3.1 For any continuous function p mapping a compact and convex set U of an 
Euclidean space into itself there is a point f such that = f. 

Note, the uniqueness of the fixed point f is not guaranteed under the conditions of 
Brouwer’s theorem. 


Theorem 3.2 Assume that arrays L = {kjk)^^nxp ® = i^iok)mxnxp consist of non¬ 
negative values and satisfy the following conditions 

(i) If >0,ll> 0,lf > 0; 

(a) Cijk > 0; 1 = 1,..., to; j = l,...,n; k = l,...,p. 

Then the system of difference equations defined by \8.1\) has at least one solution f G 

j^m+n+p 


Proof. According to the usual notation i?”*, i?", and are Euclidean 

spaces. Let us introduce vectors 

f := (x; y; z) := (a:i, ... , yi,... , y„; Zi, ... , Zp) G i?™+"+P 

and the function (p : R'^+'^+p jpn-i-n-i-p following formulae 

n p 

X] S cijuyjZk 

T^{i) = -, i = l,...,m, (3.2) 

^ E E ^ijkVj 


(f) 


m p 

iV '^2/ '^2 ^ilk^i^k 
i^lk^l _ 

jV m p 

^2 ^ijk^i^k 
i—1 k—1 


i = 1, 


n, 


E E CijiXiyj 

Tm+n+k{f) = -, k=l, ...,p. 

E E CijkXiyj 

i=lj=l 
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According to the assumptions of the theorem the denominator of the rational functions 
ifii, (fim+j and (fim+n+k are never equal to zero, hence, the vector function ip is continuous. 
Let us denote 


:= min {cijk : j = 1, n; k = 1, ... ,p}, 

Mil; = max{eijk ■ j = 1, ..., n; fc = 1, ... ,p} , i = 1, ..., m; 

fji^ := min{eijk : i = 1, ..., m; fc = 1, ... ,p} , 

Mj = max{eijk : i = 1, ..., m; fc = 1, ... ,p} , j = 1, ..., n; 

pI := min {e^jk : i = 1, ..., m; j = 1, ..., n} , 

M! = max {eijk '. i = 1, ..., m; j = 1, ..., n} , k = 1, ..., p. 


Obviously, the condition (ii) is equivalent to 


/i“ > 0,i = > 0,j = pi > 0,k = 1, ...,p. 


It also holds 


Ifpf 


<Pi< 


If Ml 


IfMf - 
l]pl 
IIMJ 

W 


^ Pm+j ^ 


, < 




IfMf 


^ V^m+n+fc ^ 


ifp 

llMl 




* = 1, ..., m, 
j = 1, ..., n, 
,k = l, ...,p, 


(3.3) 


i.e. the vector {tpi, ... prn,P 7 n+i, ... ,Pm+n,Pm+n+i, ..., v^m+n+p} bclougs to the closed 
box U C defined by inequalities (13.3|) . Therefore, the function p defined by for¬ 

mulas (13.21) maps any initial guess f° := (x°;y°;x°) into U and maps box U given by 
(|3.3I) into itself. 

Closed and bounded box U is compact and convex set in the Euclidean space and the 
function p is continuous. Hence, all conditions of the Lemma 13.11 (Brouwer’s fix-point 
theorem) are satisfied. Therefore, there exists a fixed point f of the function p defined 
by (13.21) . i.e. it holds (p(f) = f, which proves the theorem. 


Note here that seemingly restrictive conditions (ii) cannot be loosened. If, for exam¬ 
ple, Bijk = 0, then pf = p^. = pI = 0. It means that the box U is not bounded and 
the function (/^(f) is not necessarily continuous. Hence, the assumptions of Brouwer’s fix- 
point theorem would not be fulfilled. However, the condition (ii) is not a restriction in a 
real life situations. Namely, it makes sense to do the adjustments of factors by formulas 
(|2.8p only if we have statistically credible experience of all slices, which requires large 
volumes of insurance in play i.e. large exposures to risks. That is highly correlated with 
the condition (ii). If it happens anyway that for some cells it holds Cijk = 0, then rating 
actuaries need to merge neighboring slices and that way eliminate zero-risk cells. In fact, 
the slices that do not have credible experience are irrelevant in the rating process and 
actuaries merge them with neighboring slices to create slices with credible experience. 

Note also that from practical (actuarial) point of view, we need a unique solution of 
the iterative process iU), i.e. a unique solution f = (x;y;z) of the equation (p(f) = f, 
while we only proved that under conditions of the theorem there exists at least one fixed 
point of the function p. 


3.3 Uniqueness of solution: In order to discuss sufficient condition for convergence 
of the algorithm dSH]) recall the following well known fixed-point theorem, see e.g. [6] 
Section 4.1. 
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Lemma 3.3 If (p : U —>■ R‘^ is a continuously differentiable function in a convex set 
U C and if there exists a constant positive number p < 1 such that for any norm in 
it holds 

||J^(f)|| <P, Vf GU, 

then (fi has a unique fixed point f m U and the iteration U.3\) converges to f € for 
any initial guess f° chosen in U. 


In order to find conditions that will guarantee uniqueness of the solution of the func¬ 
tion (fi we need hrst to derive J^p. We again have 

f := (x; y; z) := (xi, ... , yi,... , y„; Zi, ... , Zp) G i?™+"+P. 

In order to see how the Jacobi matrix 

d{xi, . . . , Xm, yiy ■ ■ ■ y Vmy Zl, . . . , Zp) 
looks like we will introduce some block matrices. For example, 


d{yiy ■■■yVn) 


/o 


9^1 

dV2 


0 9^ 

9y2 


dvi \ 

9y„ 


9vm 
9yn / 


where 


% - 
yj If dy^ 


f \ 

eijkVjZk 

j=l k=l 


dVy II dy, tte^jkVjZk 

\j=ik=i ) 

n p 


^^lj_ \k=l 

dVj If 


X/ ^IjkZk ) X/ X/ ^ijkVjZk I X/ ^ijk^k ) S S ^IjkVj^k 


j=l k—1 


/c=l 


n p 




' n p ' 

^ijkUj^k 
k^i j 


Obviously, 

m) ^{yPm+ly ■ • ■ y V^m-t-n) ^{y^m+n+ly • ■ • y ^m-t-n-t-p) 

d{xi,...,Xm) d{yi,...,yn) d{zi,...,Zp) 

Then, the matrix represented in terms of block matrices is 

/ 0 


= 0 . 


d'p - 


V 


_ . 9 (zi,...,zp) 

0(cp^+l,...,cp^+n) PI 9(<Pm + l,...,ym + r.) 

d{x-l,...,Xm) d{zi,...,Zp) 

d(ipm + n. + l,...,<Pm + n.+p) 3(y m+n + 1, ■ ■ ■ m+n+p) 

d(xi,...,Xm) d(yi,...,yp) 


3(yi,...,cpm.) 

d{yi,...,yp) 

0 


0 


(3.4) 


(3.5) 


3.4 Sufficient conditions for convergence: To find sufficient conditions for con¬ 
vergence of the algorithm defined by (13.11) we will first work with || 

From definition (13.3|) of the box U, if > 0, by maximizing numerator and mini¬ 
mizing denominator we get 


dyj - If' 


n p 


Mf Y) zk] { Mf Y Vj zk] - [ pf Y zk] [pf Y vj E ^k 


j=l fc=l 


fc=l 


n p 




' n P 

pf E Vj E zk 

y k^l J 
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I- (Mf Mf - 

IX n 

' (Mf)" E % 


The same estimate we get in case < 0. Thus 


d'fi 

dy, 


If (Mf Mf - fiffif) 

]x n 

^ E % 


Then, from the second relation in (13.31) we get 


< 




T.V, E S- 




i=i 


i=i 


m^E M 
- 


and finally 




dy. 


< 


If (Mf Mf - Ai>f) If 


IX 


(A^f) 


n ly 

ijV V 
PI 


i = 1, ..., to; j = 1, ..,n. 


By the same token 


dtp. 


m+j 


dxi 

^^m+j 


< 


ly {MfMy - pipy) if 


n 




IX 

A^T E a7^ 

i=l 


dzk 


< 


ly {MfMy - pIpI) If 




, j = 1, ..., n; i = 1,... ,TO, 


j = 1, ...,n-,k = l,...,p 


yf E 


fc=i 


fc 

m; 


etc. Then, for every vector / = {xi, ... Xm, J/i, • ■ •, 2/n, zi, ..., 0 p} S U C 
have 


V ^ ^ ^ It {MfMf - pfpf ) _ ll 


(A^f)' 


n ,y 

..y 

PI MV 

j=l ^ 


= (n- 1) 


^ ^^lf{MfMf-pfpf) If , 

Z ^ -rTT2-= b-i) 


If {MfMf - pfpf) _ ll_ 

If {pff 

If {MfMf - pfpf) I 


n ,y ’ 
V Y' -U- 

ri M“ 

j=i J 


k=l 




A^f E 


fc=l 


k 

Mf 






A^f E 


fe=i 


fc 

m; 


etc. Thus 


< 


Af (MfMf - /r>f) 

(/if)" 


E 

/=! 

/ 




9to 


(n-1) 


E 

/c=l 


ll 


dPi 


dzk 


< 


V 


V -X 

A^i 2^ mV 

i=i ^ 


+ (p - 1) ■ 


Aif E 


fc=i 


Mi 


= '■ Po 


J 


Note, we here introduced the number pf^. Similarly we introduce Pc, 


E 


c/p. 


m+j 


dxi 


E 




m+j 


dZk 


< 
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zj (Mf M/ - 


( 


II 


{A 


(to — 1) ■ 


V 


Ml S M® 


+ (p -1) 


n 


=: Pio 


E 


^V^m+n+fc 




E 

i=i 


^^m+n+k 


dVi 


Ml S MJ / 

fe=i / 

< 


II (MfM| - AifAil) 
(M^)^ 

We also introduce 
From (dU and (1^ we get: 


/ 


- 1)- 


V 


Ml J2 Afi» 


Mf E 


=: Po 


Mj 

=1 ^ / 


Moo := max{p(„,/4,p^} . 
* 


IIEIL- 


max 


la P o m 

^ dyj azfc 




m+j 


dxi 


k=l 


dy}. 


dz^- ’ ^ 


i=l 


^V^m+n+fc 


(9x, 


E 

i=i 


^^m+n+k 


dVi 


Now, it obviously holds || < Poc- If Moo < 1, according to Lemma 1331 the function 

defined by (13.21) has a unique fixed point and algorithm (13.1|) converges to that point. 

That proves the following theorem that gives us sufficient conditions for convergence 
of the iterative algorithm given by formula m- 

Theorem 3.4 If conditions of Theorem \3.2\ are satisfied and if Poo < 1, then the function 
defined by formulas li3.S\} has a unique fixed point, and algorithm i3. 1\) converges towards 
that point for any initial guess f® chosen in . 

3.5 Criteria obtained by : Similar results can be obtained by means of 

II We start from (II. 4p . Then 


IE 111 “ 


E 

i=i 
It holds 


max 

i,j,k 


dip, 


m+j 


dxt 


k=l 


E ^'Pm+n+k iXPi 

r)^: ’ f)'}l . • ^ 


dXi 




E 

1=1 


dp. 


m+j 


dxi 


dy, 


E 

/c=l 


fc=l 


^I^m+n+k 


^I^m+n+k 


dx. 


dy, 


< 


;E 


9m* 


dzk 


E 

i=i 


dp. 


m+j 


dzk 


< 


I_[ ^ II {MfMy -pipy) ^ II (MfM| -Pipi) \ 


Mf E 


Mf 


E 

. 1=1 




E 


{Plf I 


, TO. 


Note, we denoted by p\ the right hand side of the above inequality. Similarly, we 
define 


1 L 

Ml := - - 


If {MfMr - EMf) I Ik mMl -plpl) ] ^ _ 


1=1 

k If 
P i ■■= 


_ ) q yiviiivi^ Pi Pi) 

V ^ I ^ ^ 

Ml E H 


tflf (Mff ' [plf 


If {MfMf - pfpf) ^ 1] {MfMy - p{py) 


fjL fcr'f jzpi (+tf 


fc = 1,... ,p 


Ml ^ M, 

fc=l ' 
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and 


Pi ■= . 

As II ^ Pi, the following statements follow; 

Theorem 3.5 Assume that the conditions of Theorem \S.2\ are satisfied. 

(i) If Pi < 1, then the function defined by formulae 2\} has unique fixed point, and 
algorithm i3.1\) converges towards that point for any initial guess chosen in D°. 

(a) If p := min {pi, Poo } < 1, then the function (p given by \3.2\) has a unique fixed 
point, and the algorithm a converges towards that point for any initial guess 
chosen in D^. 


3.6 Simplified criteria for uniqueness; If we now introduce 


p ;= min{eijk : i = 1,... ,m; j = 1,... ,n; k 
M := max {cijk '■ i = 1,... ,m;j = 1,..., n; fc 


then we can estimate 


1 1 

- < - 

n — n iv ,,y 

fe=i j=i ^ ^ 


Km 
< A — 
Lp. 


1 ,... ,p} 


and 






< 


If M (m2 - IV 




pfiL 


i = 1, 


to; j = 1, ..,n 


etc., then all derivations and statements simplify significantly. Hence, repeating the proof 
of the Theorem 13.41 gives; 


Corollary 3.6 A ssume that the conditions of Theorem \3.2\ are satisfied and 

'^oo ■ L 

1 . 

Then the function defined by formulas i3.2\) has a unique fixed point, and algorithm i3.1\} 
converges towards that point for any initial guess f chosen in . 


' (n-i)qi\ 
If 


(p-i)ifd. , (p-i)ifil , 


{m—1)1^1^ 

If 


{n-l)liq 




< 


Repeating the proof of the Theorem 13.51 with M and p we get the following very 
simple criteria; 


Corollary 3.7 Assume that the conditions of Theorem \3.2\ are satisfied and 


M (m 2 - p2) 

ri ;= -5- -uidx 

p3 


II 


If II 

ir If 


1. 

If’ 



< 1 . 


Then the function defined by formulas i3.2l) has a unique fixed point, and the algorithm 
converges towards that point for any initial guess chosen in . 

If r = min {ri, Too } < Ij then the function ip given by i3.0i) has a unique fixed point, 
and the algorithm mD converges towards that point for any initial guess f® chosen in 
D°. 


Obviously, poo < ^oo and pi < ri. Therefore, the simple criteria of convergence given 
by corollaries are weaker than the criteria given by theorems 13.41 and 13.51 
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3.7 Note on generalization: If we have N>3 risk classification variables class, 
territory, ..., industry, i.e. industry is rather than third classification variable, 
then we have for class i 

n p 

Zf ^ i = l,2,...m, (3.6) 

j=i fe=i 
n p 

'.= 'y ( ... y ( Cijkyj ■ ■ ■ Zk, 1 = 1 , 2 ,..., 771 , 
j=l k=l 

meaning that only summation by index i, corresponding to factor x, would be omitted. 
Similarly, for all other risk factors y = (jji, ..., y„), ..., z = (zi,..., Zp), summations 
by the corresponding indexes are omitted in the formulas of the form (|3.6I) and (13.71) . 
Then, all other formulas, including (j2.4|) . (12.81) . (j2.9l) would be generalized accordingly. 
For example, the function (13.21) will be defined now by 

n p 

,2. S • • • S 

l~ j=l k=l 

lx n p 

^ E • ■ • E eijkVj 

j=i k=i 

and similarly ..., ipm+n+...+k- Then generalizations of all the statements will be 

cumbersome, but straightforward. 


...Zk 

-, 1 = 1, . . . ,TO 

■ - -Zk 


4 Application to Multi-species Leslie-Gower Model 

4.1 Leslie-Gower Model: Competition of d > 2 species in an ecological system is 
modeled by the following Leslie-Gower system of difference equations 

bixi 

xi = - - 

1 -I- CiiXi -I- ... -I- CidXd 


Xd 


_ bdXd _ 

1 + CdlXi -I- ... -I- CddXd 


(4.1) 


Hence x = (xi, ..., X2) again denotes the next iteration of x = (xi, ..., X2). It is natural 
to assume 


bi ^ 0, Cii ^ 0 , Cij ^ 0 , 


h J 


= 1 , 


d, (xi, ..., Xd) € := [0, 00 ) X .. .X [0, 00 ) C R'^ 


Here Cij = 0, i 7 ^ j means that the species i is not affected by the species j. The 
models for two species (variables) are extensively studied. Study of systems with d > 3 
species is more complicated for many reasons. For example, some research for two species 
reduces at the end to a quadratic equation. Therefore, the same methodology cannot be 
extended to higher order systems, because we would end up with algebraic equations of 
higher order that cannot be solved. 

The corresponding Beverton-Holt function ip : R‘^ R'^ can be written in the form 


(pi{x) 


6 ,; a 


1 “t“ CiiXi ^ ] CijXj 


,i = l, 


d. 


(4.2) 


4.2 Necessary conditions for convergence of Leslie-Gower model: 
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Lemma 4.1 (i) If system o converges to a solution 

X € := {x € R‘^ : Xi > 0, yi = 1, ..., d} , 

then 


> 1, V* = 1, ..., d, 


and 


Range 


(ii) The matrix 


Cll 


Cdl 


Cid {bi - 1) 
Cdd {bd — 1) 


= Range 


Cll 


Cdl 


Cld 


Cdd 


C = (c„) 


'^3 3 dxd 


(4.3) 


(4.4) 


(4.5) 


is invertible if and only if the fixed point of p, which is also a solution of the system 
is uniquely determined. In that case J^.dl ) and o also hold. 


(in) If Cij = 0, \/i j, then x reduces to 


K := 


bi-1 bd -I 
Cll ’ ’ Cdd 


e 


(4.6) 


Proof For every solution x := (si, ... , Xd) cD^ of the system (14.11) . it follows 


Ciia^i + ... + CidXd = bi - I, i = I, ..., d. (4.7) 

Now the first statement in (i) follows from the positivity of the left sides of (lOl . The 
second statement in (i), as well as the statement (ii), follows from the well known facts 
of linear algebra. The remaining statements are obvious. 


The two necessary conditions for convergence of the LG system (El) are listed in the 
statement (i) of Lemma l4.11 and necessary and sufficient condition for the uniqueness of 
the solution is given by (ii). As we know from linear algebra, equation (14.4p is necessary 
and sufficient condition for existence of the solution of the system (ED- 

Remark 4.2 A multi-species Leslie-Gower models were studied in 0/ and The fol¬ 
lowing model was proposed in 

^ t^iR-iXi . , 

Xi = „ , . -TT--, 1 = 1, ..., d 

Ri \t^i IjXi -\- / . CijXj 

jAi 

where Hi > 1, Vi = 1, ..., d was assumed. The weak inter-specific competition was 
characterized by the condition stated there as: ’’Coupling terms Cij, i j are sufficiently 
small. ” 

It is easy to prove that relations between coefficients bi and Cij in o and the coef¬ 
ficients pLi, Cij and Ki in JBl ore the following 

Pi — bi, Ki — ^ , Cij — KiCij, i j, i,j — 1, . .., d 

The notation in m is as ours except that authors there use different letters and they 
set variables to obtain the diagonal coefficients equal to 1; in our notation that would be 
Cii — 1 . 
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Note that in Lemma |4.1l we did not assume weak inter-specific competition and we did 
not assume 5^ > 1, Vz = 1, ..., d; we obtained that as a necessary condition for con¬ 
vergence of any system of the form (HU. We will later deal with weak inter-specific 
competition, see our assumption (I4.8|l . 


4.3 Existence and uniqueness of a solution of a multi-species LG model: 


Let us briefly distinguish the results obtained in our Theorem 14.31 from the corre¬ 
sponding results in [8] and [7]. 

The assumptions in [5] for the existence and uniqueness of the solutions of the system 
(|4.1I) were established by relations (13), (14) and (19). Our assumptions in Theorem 14.31 
are simpler and more general and our proofs are simpler as well. For example we do not 
need to use dynamic reduction technique. Regarding [7], one of our assumptions, it is 
(j4.8l) below, is equivalent to (24) in [7]. 

The real difference in the following Theorem 14.31 makes assumption that the matrix 
C = is invertible. That condition is the most general condition for uniqueness 

because it is also a necessary condition for uniqueness, c.f. Lemma |4.II fiil. In addition, 
that assumption enables the following surprisingly short proof of the theorem. 


Theorem 4.3 Let the function ip : ^ defined by equations satisfy condi¬ 

tion |/).d[) and 


Cij — <bi-l,'ii = l, 


For any 


we can introduce the box 


hi < 


-EC: 


^3 
V Cjj 


d. 


(4.8) 


(4.9) 


U := j {xi, ...,Xd) & R+ ■■ hi < Xi < —,\/i = 1, 2, ..., 

Cii 

(i) Then the funetion ip has at least one fixed point x := (xi, ... , 

Every such point is also a solution of the system 

(ii) If additionally matrix C = icij)^^^ is invertible, then the solution is uniquely de¬ 
termined, and it is a solution of the linear algebraic system O- 

(Hi) The algorithm 0]7]) converges to that solution independently of the selected initial 
condition x^ = (xfj ■ • ■ 1 3:°) ^ ■ 


Xd) in the box U. 


Note, here we assumed weak inter-specific competition by condition (|4.8I) . 


Proof, (i) In the first step we will prove that the function p maps the box U into 
itself. 

For every x € it holds 


Pi (x) < 


biXi 

CiiXi 



(4.10) 


Therefore, it remains to prove 

hi < pi{x) ,yi = 1, ..., d. 
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Assume x € U. From the definitions of U it follows: 


(Pr (x) := 


biX, 


1 “t“ CiiXi F y ' 


> 


bihi 


1 F Ciih'i F y CijXj 




bihi 

> - — - j— > 

1 F Ciibi F X/ ’'ij pjj 


^iihi F ^ij f,.. 

Only the last inequality is not obvious. However, it easily follows from the assumption 

Hence, the continuous vector function p maps the box U, a convex and compact set 
in ^ into itself. As before, by means of Brouwer’s Fixed Point Theorem we conclude 
that the continuous function p has at least one fixed point in U. 

(ii) Assume additionally that the matrix C = is invertible. Then any solution 

X := (xi, ... , Xd) G U of the system (14.111 . which exists according to (i), satisfies 
all conditions of the Lemma l4.II (iil. Therefore, the solution is uniquely determined. 

(in) According to (I4.10|l . already the first iteration satisfies 

Xi < —, Vz = 1, ..., d. 

C-ii 

Then we can simply select hi in the definition of the box U sufficiently small to satisfy 
both, (14.9p and /li < Xi, Vz = 1, ..., d. Then the statement (iii) follows from previous 
statements of the theorem applied to that particular box U. 


This theorem practically assigns a corresponding linear algebraic system of the form 
S3 to each system of d difference equations of the Leslie-Gower’s type with weak inter¬ 
specific competition. That way the study of the system of difference equations of the 
Leslie-Gower’s type is simplified to study of the corresponding linear algebraic system. 
More formally the above results can be summarized in the following characterization. 

Corollary 4.4 Assume that the system ( S3 satisfies condition 

Then there exist a unique solution x S U C if and only if the following two 
conditions are satisfied: 

(a) > 1, Vz = 1, ..., d. 

(b) The matrix C = (cij)^^^ is invertible. 

That solution is independent of the initial guess x°. It can easily be calculated as a 
solution of the linear regular algebraic system ra and it is also the unique fixed point 
of the corresponding Beverton-Holt function p. 
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